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GOALS OF THRUSTER MODELLING 

There are currently no thruster modelling codes that can be integrated with 
P»wer system codes for full propulsion system modelling. SI eriiSSttS 
Wri,le " from » "*•*»« a'«"e" viewpoint. assuX tL uTr u 
m,!H ii™ *ir an f - y “ S J °n ,hrus ‘er performance alone. The goal of*the present 
modelhng effort is to develop thruster codes that model performance and scaling as 

phy^cXValter and ^ infUtS ' ra ' her ,h * n » o' more elemental 

level parameters of interest are performance, such as snerifir i mniI i 0fl 
and efficiency; terminal characteristics, such as v”^ge or cu^em anZ^ 

,n *P uls * an< J efficiency couple with mission analyses, while terminal 

Hfp#im ter ? CS a ?- W nte8ration with Power systems. Additional information on 
lifetime and operating may be required for detailed designs. 
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APPROACH TO MODELLING 

• Two Thiuster Types 

- Inert Gallon 
-Hydrogen MPD 

• Focus on macroscopic characteristics, not microscopic 

• Use fundamental physics (where possible), empirical 
data 

• Benchmark results with existing experimental, design 
results 

• Currant results do not represent optimized designs 


APPROACH TO MODELLING 


For this initial effort, the two thruster types with the strongest development 
background are being modelled: the Magnetoplasmadynamic (MPD) and Ion 
Thrusters. The emphasis is on modelling these devices as systems; that is, to focus on 
the macroscopic system level parameters such as power, thrust, specific impulse , 
rather than on the microscopic parameters such as electron temperature, ionization 
fraction, and plasma instabilities. Where possible, the fundamental physics of the 
concept are used, to provide as close an understanding of the underlying processes as 
possible. Where understanding is incomplete, or too complex for productive system 
analysis, empirical results have been used. For example, applied field MPD thruster 
thrust generation is based on experimental measurements, rather than an analytical 
model. 

As these models are developed, they are and will be compared to experimental 
data and point studies. 
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MPD Thruster 

i «r»«J h K thruster accelerates a plasma propellant through the electromagnetic 
^ . ! z ^ ody forc ?* The system considered in this modelling activity is a cylindrical 

coaxial thruster, with an external anode and central cathode. Acceleration k n?oS 
through the interaction of radial and azimuthal currents with both the self-induced 
(azimuthal) and applied (axial and radial) magnetic fields. The annlied fi*M ■« 
generated by a solenoidal coil located externally of the anode. The majority of the 
thruster's waste heat has been observed to be deposited in the anode, reauirinl a 
radiator to reject this energy to space. In this design, the radiator is a se^ Juthium 
heat pipes conductivity coupled to the anode and transferring the heat from the anode 
surface to a surrounding circular graphite surface. 

. .... Constraints on MPD thruster operation are cathode lifetime due to mass loss, the 
ability to reject the anode heat, material temperature limits, and the cooling 5 the 
hyperconducting magnet coil, which operates at 21 K ’ a coo,,n & or the 

Inputs range from performance requirements to some system design parameters. 
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Performance: 2.5 MWe» 5000 s,t| 

■ 0.5, BO « 0.5 T 




Design* 

Model 



Anode Radius (cm) 

IS 

15 



Cathode Badka (cm) 

2.5 

2.6 



Anode Length (cm) 

30 

30 



Cathode Length (cm) 

10 

10 



Current (kA) 

10 

8.5 



Voltage (V) 

250 

295 



Magnet Current (A) 

2300 

2471 



Anode Temperature (K) 

1400 

1861 



Anocto FaH Volta* (V) 

25* 

90** 



Radtator Area (m2) 

1.1 

4.4 



Maas (kg) 

-337* 

993 



Maas wfo Radiator (kg) 

-132* 

117 



♦Myers, et *L "Mullimcgawatt MPD Thruster Design 
Symposium on Space Nuclear Power Systems, 1992 

Considerations," 9th 



1 * Assumed **Bstimated 

ffPnxn a Related Study of a Flared Anode Thmster J 



MPD Thruster Model Benchmark 


An initial benchmarking of the code in terms of system level parameters has been 
performed. The point design is actually a combination of results from two references: 
"Multimegawatt Electric Propulsion System Design Considerations," AIA A 90-2552; 
and "Multimegawatt MPD Thruster Design Considerations," in the 9th Symposium on 
Space Nuclear Power Systems, January, 1992. MPD thruster mass was taken from the 
first reference, which actually used a flared anode, with an inital anode radius of 15 cm 
flaring to 30 cm at the exit. The second reference is a cylindrical anode of 15 cm radius. 
The second reference was used for input data to the MPD model. 

In terms of terminal characteristics and magnet design, the model results are 
reasonably close to the point design. Such differences that do exist are due to differences 
in assumptions of applied field thruster performance, and could be remedied through 
better empirical parameters in the model. 

Model results differ primarily in terms of radiator mass. This is because of the 
difference in anode heating between the two cases. The reference case assumed a low 
(25 V) anode drop, whereas the MPD model estimates a 90 V drop. This difference 
shows up in both the radiator size and the anode temperature. An improved model of 
MPD thruster loss mechanisms will be required to resolve this difference. 
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MPD Thruster Model 
Effect of Magnetic Field on Thruster 

Design 

The impact of the applied field upon thruster design is shown in this figure. 
Increasing the applied field increases its contribution to accelerating the propellant, 
reducing the need for the self field thrust component. This results in a decrease in 
anode radius, for conditions of constant power and efficiency. This effect is seen to 
become less marked at higher fields, indicating that there may be maximal field 
strength for MPD thruster operation. 

This result indicates one benefit of the model: previously, scaling of the thruster 
with field strength had not been addressed on a parametric basis. Instead, a single 
design point of field strength and anode radius was selected. It should be noted that 
this anode radius is also consistent with anode heat rejection and heat conduction 
constraints. 
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MPD Thruster Model 
Effect of Magnetic Field on Thruster 

Mass 


fi*M Jbe scaling of thruster system (anode, cathode, magnet, radiator) with annlied 
field is shown here. The result indicates a region of field strenShs wST JEli 
raster system mass. In the present model, radiator mass is a dominant segment of the 
m,n 7 m ®ass Point is due to a trade off in decreased anode and maraet 
0SSCS at hi8her ne,ds * This behavior is dependent uSmthe 
f Um P t, ® ns ’ curren *ly an area of experimental and theoreticalinvesUgation 

amranhi* * anode n,odel W,H ensure the minimum mass point. The MPD model is 
amenable to incorporating such changes as they become necessary. 
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MPD Thruster Model 
Effect of Efficiency on Thruster Mass 

The dominance of radiator mass in the overall system mass is seen in this 
calculation of thruster mass for varying efficiencies. Increased efficiency is simply 
decreasing the amount of waste heat delivered to the anode. Additional effects due to 
thruster or magnet radius are subsumed in the radiator effects. 
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Ion Thruster 

p.«j^r, s v h n^rs at zzsrrtt?, r/j 

downstream from a discharge chamber in which the plasma is generated 8 Proliant. 

"*.?V TOr | 8ases k W»n, and «rgo“proS„ t rtoi«dC„t 

upon the specific impulse and efficiency required. aepends 

. th " ,sters operate under several constraints. The primary limit is the snace 

constrain/* 1 U * ,0n ,on beam density. In addition, numerous engineering level 
constraints upon power density exist, such as grid lifetimes. These considerations »r*» 

pHH r?Ir of P rope ‘ ,ant and operating conditions. Of the constraints listed here all but 
grid lifetime have been addressed in the thruster model to date. 

ratio S This 52 L b T d ?“ en ? ineerIn 8 concerns, such as the span-to-gap 

„? n tl l f , the thruster g nd length (the span) to the inter-grid spacing 

this rSJ^rTXtati^' de,0rma “ on - ,hcre b « P™*** -PI™ limit to 
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Inn Thruster Example 


Inputs: 

Propalant 


top 

Max Span-to-Qap 
Max Voltage 


Xa, Kr, A 

100 kWa 

3000 - 10000 $ 

880 

3000 


Outputs: 

No. Thruntert 
Beam, Dtecharga 
Current, Voltago 
EHJcwky 
Area 

Thruster Mats (TBD) 
Lifetime (TBD) 


Ion Thruster Example 


A sample case of a 100 kWe ion propulsion system has been assessed for this 
presentation. Inputs are shown above. The ion thruster model was used to calculate 
system parameters and operating conditions that both met the input requirements and 
satisfied the constraints. The thruster model will ultimately calculate thruster masses, 
as does the MPD model. 
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Ion Thruster Model 
Thruster Performance 


nrm JtZ.[ hr 'i S l er p * rfol ' mance (efflcency, specific Impulse) is shown for all three 
Jj-iV 8 These results are comparable to experimental data for 30 or SO cm 
diameter thrusters operated at Lewis Research Center. It -i-v ld be noted fhm III 

outpuTonh^°modeT ated ^ ' hrUS ' er dimensio,,s; n, * her - ,hrus «* «•«■>* was an 
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Ion Thruster Model 
Scaling of Thruster Size with Isp 

Thruster scaling is shown for the three propellants. Total grid area is the area 
required to process 100 kWe of power, although the number of thrusters changes with 
specific impulse. The model predicts greater power densities at higher specific impulse, 
as is seen in experiment. The behavior of these data may change after grid lifetime 
constraints are imposed. 
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Ion Thruster Model 
Relation of Isp to Span to Gap Ratio 

. .™ e re fl u,ret | span-to-gap ratios for each operating point are shown As nower 
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Summary of Progress 


• Development of system level thruster models is underway 

• MPD and Ion thruster models are aimed at Integration with 
power system and trajectory codes for trade analyses 

• An initial ”End-to-End H MPD thruster model has been 
developed and compared to point designs and experiment 

• An initial ion thruster model is under development 

• Both models have the potential for additional refinement 
through the use of more complex computer codes where 
needed; however, simplicity to desirable for rapid trade 
studies 


Summary of Progress 


This presentation is intended as a status report on thruster system modelling 
efforts currently underway at Lewis Research Center. An evolutionary approach is 
being taken in developing these models. Refinement of the codes and their component 
subroutines is expected in the coming months. First order modelling has provided some 
initial insights into thruster behavior and requirements for effective implementation. 
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Further Work 


J n add,t,on t0 completing the ion thruster code lifetime and mass models, several 
lmp ™ ve ™ ent of . both codes are evident. The impact of the MPD power toss 
models upon thruster design emphasizes the need for a better understanding either 
theoretical or empirical, of dissipation in the MPD thruster. Furtter refinement of thl 
radiator model is required for effective system design. refinement of the 

In both codes, the potential for internal optimization of certain thm^r 

i0B ,hr “ ster « rid «—• - «* 

. , Pe "J. 08 * im P° rtant at this Stage is that thruster system models are beimr 

prwe^f involved^ ^ ana,ySlS ^ providin 8 some understanding of the physica? 
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